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> ' ABSTRACT 

r- : 

Aims. A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is employed to investigate the properties of 
SNR RX J1713. 7-3946. 

Methods. Observations of the nonthermal radio and X-ray emission spectra as well as the H.E.S.S. measurements of the very high energy 
' y-ray emission are used to constrain the astronomical and the particle acceleration parameters of the system. 
, Results. Under the assumptions that RX J17 13.7-3946 was a core collapse supernova (SN) of type Il/Ib with a massive progenitor, has an age 
■ of ~ 1600 yr and is at a distance of a 1 kpc, the theory gives indeed a consistent description for all the existing observational data. Specifically 
it is shown that an efficient production of nuclear CRs, leading to strong shock modification, and a large downstream magnetic field strength 
Sd ~ 100 fiG can reproduce in detail the observed synchrotron emission from radio to X-ray frequencies together with the y-ray spectral 
Q characteristics as observed by the H.E.S.S. telescopes. Small-scale filamentary structures observed in nonthermal X-rays provide empirical 
j confirmation for the field amplification scenario which leads to a strong depression of the inverse Compton and Bremsstrahlung fluxes. Going 
CZ2 . beyond that and using a semi-empirical relation for young SNRs between the resulting CR pressure and the amplified magnetic field energy 
upstream of the outer SN shock as well as a moderate upper bound for the mechanical explosion energy, it is possible to also demonstrate the 
^> actual need for a considerable shock modification in RX J 1 7 13. 7-3946 . It is consistent with RX J1713. 7-3946 being an efficient source of 
nuclear cosmic rays. 
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1. Introduction 

RX J17 13-3946 is a shell-type supernova remnant (SNR), lo- 
cated in t he Galactic plane, that was discove red in X-rays with 
ROSAT (Pfe ffermann& Aschenbachlll996t). Further stud y of 
this SNR with the ASCA s atellite bv lKovama et alJ fl997) and 
later bv lSlane et alJ i 19991) have shown that the observable X- 
ray emission is entirely non-thermal, a nd this property was 
confi rmed in later XMM observations ( Cassa n>Chenai et alJ 
2004). 

Given these characteristics it is somewhat surprising that 
the radio emission is weak. In fact, only part of the shell could 
be detected in radio synchr otron emission up to now, with a 
poorly known spectral form jLazendic et all2004T) . 

RX J1713-3946 was also detect ed in very high en- 
ergy y-rays with th e CANGAROO jMuraishi et alJ [2000: 
Eno moto et all 120021) and H.E.S.S. jAharonian etaD 12004 
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2005) telescopes. Especially the latter observations show a 
clear shell structure at TeV energies which correlates well with 
the ASCA contours. 

This intriguing situation has prompted us to theoreti- 
cally investigate the acceleration of both electrons and pro- 
tons in detail, using a nonlinear kinetic theory. It couples the 
particle accelerati on process with the hyd rodynamic s of th e 
thermal gas JBerezhko etaT]ll996t iBerezhko & Volkll2000al) . 
However, in comparison with the other SNRs (SN 1006, Cas A, 
and Tycho's SNR) that were successfully described within 
the framework of this the ory and finally allowed a predic- 
tion of the T-ray flux JBerezhko et alJl2002l: IVolk et alJl2002t 
I Berezhko etafll2003albt jBerezhko & Volklbo04al IVolk et alJ 
120051: iKsenofontov et aljbo05l) . the case of RX J17 13.7-3946 
puts considerably more obstacles in theory's way. First of 
all, such decisive astronomical parameters as source distance, 
present expansion velocity and age are not well known. We 
shall use the values d — 1 kpc and t = 1612 yr (see below). 
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In addition, the lack of knowledge of the spectral shape of the 
radio emission makes it difficult to derive - from synchrotron 
spectral observations - two determining physical quantities: the 
effective strength of the magnetic field and the proton injection 
rate into the diffusive acceleration process. This technical dif- 
ficulty adds to the uncertainty about the type of progenitor star 
and its possible influence on the circumstellar medium. At first 
sight it appears therefore impossible to theoretically predict the 
TeV y-ray emission. We shall argue nevertheless that the ob- 
served overall synchrotron spectral shape, from radio frequen- 
cies to the X-ray cutoff, and the small-scale filamentary struc- 
tures in the nonthermal X-ray emission of RX J17 13.7-3946 
are consistent with efficient CR acceleration associated with a 
considerably amplified magnetic field. These properties allow 
in addition a consistent fit of the observed TeV y-ray spectrum. 
It is strongly dominated by 7r°-decay emission. 

Assuming that all young SNRs obey a semi-empirical lin- 
ear correlation between the amplified magnetic field energy 
density and the pressure of accelerated CRs, and that the me- 
chanical explosion energy of RX J17 13.7-3946 has an upper 
bound of 2.5 x 10 51 erg, also makes it possible to approximately 
determine the proton injection rate and thus to predict the time 
evolution of the CR energy content of the SNR and the nonlin- 
ear backreaction of the CRs on their own acceleration. 

The existence of a compact central object, the X-ray point 
source 1WGA J 17 13.4-3949 - most probably a neutron star 
- and the properties of the observed nonthermal emission of 
RX J1713. 7-3946 strongly indicate that the source is a core 
collapse supernova (SN) of type Il/Ib which exploded into the 
wind-blown bubbl e from a massive progen itor star. This was 
also argued bv lCassam-Chenai et alJ d2004l) . The main fraction 
of the remnant volume contains diluted hot bubble gas, whereas 
at the current evolutionary phase the outer SN shock already 
propagates into the dense shell of the swept-up ambient inter- 
stellar medium (ISM) gas. 

The remaining, mainly quantitative uncertainties in the de- 
tailed properties of the circumstellar environment still do not 
allow a unique theoretical solution. Therefore our primary aim 
is to construct a theoretical model that is consistent in itself and 
with all available measurements. And this is indeed possible. 
Nevertheless, we show that also the prediction of a highly effi- 
cient nuclear CR production can be cogently argued. In short, 
we demonstrate that the existing data are consistent with very 
efficient acceleration of CR nuclei at the SN shock wave which 
converts a significant fraction of the initial SNR energy con- 
tent into CR energy. We demonstrate in addition that the avail- 
able morphological evidence in hard X-rays implies that the 
leptonic channels of y-ray emission (Inverse Compton emis- 
sion as well as Nonthermal Bremsstrahlung radiation) should 
be strongly suppressed, and that the calculated hadronic y-ray 
emission is consistent with the H.E.S.S. spectrum at TeV ener- 
gies and with the EGRET upper limits in the GeV range. The 
further going prediction of efficient acceleration of CR nuclei 
on the basis of two additional semi-empirical assumptions is 
left to the Discussion section. 



2. Physical parameters of RX J1 71 3-3946 

Following the estimates o f iKovama et alJ Jl997t) from ASCA 
X-ray observations, and of lFukui et alJ d2003h from NANTEN 
CO line measurements, we adopt a distance d — 1 kpc. For 
the present angular sources size of 60' this leads to a SNR 
radius R s m 10 pc. Even though a considerably larger dis- 
tance d = 6 kpc to the SNR RX J1713.7-3946 , estimated 
in e arlier stud ies (Slaneetal. 1999), is still being d iscussed 
(e.g. ICassam-Chenai et alJl2004 iHiraga et alJl2005h . the re- 
cent study of mol ecular clouds toward SNR RX J1713. 7-3946 
(Moriguchi 2005) "strongly supports" the distance of 1 kpc. 

A rather simple theoretical consideration shows that the 
SNR should be quite young in order to reach such a large size 
and to be at the same time such a bright source of nonther- 
mal X-rays. The observed energy flux of no nthermal X-rays at 
energy e v = 4 keV is vS v ~ 100 eV/(cm 2 s) JSlane et alii 999) 
which for the distance d — 1 kpc translates into the synchrotron 
luminosity vL v = And 2 vS v s 2x 10 34 erg/s. Here v is the emis- 
sion frequency, L v and S v are the spectral luminosity and flux, 
measured at the dist ance d, re spectivel y. Accor ding to the the- 
oretical estimates of lBerezhko & VoTkl J200 4b) such a high lu- 
minosity is possible only if the shock speed is not very low, 
V s > 1500 km/s. Otherwise the electron cutoff energy becomes 
so small that it implies a very low X-ray luminosity. This leads 
to a restriction on the SNR age for an assumed uniform exter- 
nal ISM. In such a case the present epoch of SNR expansion is 
described by the Sedov solution: 

fl s = (1.2xl0 24 £ sn /MsM) 1/5 r 2/5 , (1) 

where £ sn is the SN explosion energy, Njsm = Pism/'«p is the 
interstellar medium (ISM) number density, pism is ISM density 
and m p is the proton mass. Since the shock radius and the shock 
speed are related as V s = 0.47? s /f the above condition gives a 
restriction for the SNR age t < 2.7 x 10 3 yr. If we assume the 
standard value £ sn = 10 51 erg for the SN explosion energy and 
use the observed SNR radius R s — 10 pc, Eq^will give us 
Msm < 0.3 crrT 3 . In such a diluted ISM the maximum possi- 
ble TeV y-ray, flux eFJe) = 10(AW 1 cmr 3 )(d/l kpc)" 2 eV 
cirT 2 s _1 JBerezhko & Volkl2000bl) and the expected nonther- 
mal X- ray flux vL v = 10 33 (tf ep /10- 4 ) erg/s JBerezhko & Volkl 
2004b) would be significantly below the observed values. Here 
K ep denotes the electron:proton ratio in the energetic particle 
population. 

This contradiction can be resolved if SNR RX J1713.7- 
3946 is the result of a core collapse SN of type Il/Ib which 
exploded into the adiabatic, very diluted bubble created by the 
wind of a massive progenitor star. This is also required by 
the probable existence of a compact rem nant in the center, the 
neutr on star 1WGA J1713.4-3949 (e.g. ICassam-Chenai etafl 
20041) . According to stellar wind theory (e.g. IWeaver et alJ 
ll977tlChevalierll982t) . the bubble is bounded by a dense, mas- 
sive shell which consists of the swept-up, cooled ISM ma- 
terial. Note that Isiane et alJ lll999l) . iFukui et alJ J2003) and 
ICassam-Chenai et alJ J2.004I) came to the sam e conclusion 
based on different types of arguments (see also lEllison et alJ 
2001). In a very general study the combined SNR dynamics 
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and diffusive shoc k accelera tion in wind-SNe has been investi- 
gated bv lBerezhko & VollJ fcOOObl) . 

Progenitor stars of core collapse SNe, which have intense 
winds, are massi ve main-sequ ence stars with initial masses 
Mi > 15M (e.g. lAbboilll982l) . In the mean, during their evo- 
lution in the surrounding uniform ISM of gas numb er density 
Pn = m n M.SM, they create a bubble of size (Wea ver et all 1 977: 
IChevalier&Liang ll989) 

R sh = 0.76(0.5MV£4 / po) i/5 j (2) 

where M is the mass-loss rate of the progenitor, V w is the wind 
speed, and f w is the duration of the wind period. This bubble 
could be adiabatic, and then have very low gas density, or it 
could be modified by mixing of the hot bubble gas with sur- 
rounding shell material. The large size of RX J1713. 7-3946 and 
its low age require an adiabatic bubble, which we shall adopt 
for the sequel. 

In order to determine the SN shock dynamics inside the 
shell we model the gas number density distribution N g as a 
constant in the bubble and as a power law in radius in the shell. 
This leads to the following radial profile 

N g = N b + (r/R sh ) 3( ^- [) N sh , (3) 

where A^ s h = o" sh /y ISM is the peak number density in the shell, 
is the gas number density inside the bubble, typically very 
small compared with the shell density, and <x sh = A^h/MsM is 
the shell compression ratio. We emphasize that the compres- 
sion ratio cr sh can exceed the adiabatic upper limit of 4 as a 
result of the radiative cooling in the shell. 
The mass of the bubble 

M b = (47r^ h /3)m p jV b (4) 

is rather small, Mb < M , in the case of m oderate progenitor 
masses M, < 2OM fe.g. lChevarier & Liandl 9*891) . whereas the 
shell mass 

M sh = 4nN sh m p dr^ir/R^-^ = (4KR 3 J3)Nis M m p (5) 
Jo 

amounts to a few thousand solar masses (see below). 

During SNR shock propagation through the adiabatic bub- 
ble, only a small fraction of the mechanical explosion energy 
is therefore given to gas of stellar origin. The main fraction of 
the explosion energy is deposited in the shell. 

Here we use the gas number density distribution N g (r) = 
p(r)/m p in the form 

N g = 0.008 + l.l[r/(10pc)] 12 cm" 3 (6) 

which gives a consistent fit for all existing data for SNR 
RX J 17 13.7-3946 . Such a distribution corresponds to a bub- 
ble with cr sh - 5, Nb = 0.008 crrr 3 , and a range 17.6 < 
R s h < 19.9 pc, created by the wind of a main-sequence star 
of initial mass 15M < M, < 2OM in a surrounding ISM 
of number density 110 < Njsm < 500 cm -3 , respectively 
JChevalier & Liand fl989). Such a high ISM density, required 
to produce a bubble of such a small size, is quite consistent with 
the value jVism = 300 crrr 3 derived from the analysis of the X- 
ray emission of SNR RX J1713.7-3946 JCassam-Chenai et alJ 



2004). This implies that the bubble is formed inside a molecu- 
lar cloud; a discussion of the CO emissio n from t he neighbour- 
hood of the SNR is given in the paper bv lFukui et alJ f2003). 

For the SNR age we use the value t = 1612 yr , which is 
consistent with the hypothesis of lWang et alJ dl997t) . based on 
historical Chinese records, that RX J1713-3946 is the remnant 
of the AD393 "guest star". 

We adopt the following values of SNR parameters: E sn — 
1.8 x 10 51 erg, M e j = 3.5M . As shown below, these parameter 
values lead to a good fit for the observed SNR properties. For 
a better determination of the explosion energy we would need 
the value of the shock speed V s which is unfortunately is not 
measured. 

2.1. Magnetic field strength 

Whereas the above physical parameters - explosion type, ex- 
plosion energy, progenitor mass loss, circumstellar density, age 
and distance - are part of an internally consistent picture whose 
detailed quantitative specification allows a consistent fit to the 
observed X-ray synchrotron and y-ray emission, the magnetic 
field strength plays a specific role in this discussion. The reason 
is that its strength decides the relative role of electrons and nu- 
clear particles in the y-ray emission. In the following we shall 
therefore go into the question of the field determination in some 
detail. 

As indicated in the Introduction, there are two obvious 
ways to determine the effective ma gnetic field in a SNR, a nd we 
have used both of them in the past JVolk et al.l2005HBerezhkol 
l2005t) for SN 1006, Cas A, and Tycho's SNR. 

The first method analyses the overall synchrotron spectrum. 
If, first of all, the radio continuum spectrum is softer than the 
test particle limit, i.e. with a photon index larger than a = 0.5, 
then this is attributed to the radiation of electrons that have been 
accelerated at the subshock of an overall nonlinearly modi- 
fied shock. This requires (i) an effective magnetic field strength 
which is above a certain limit and (ii) that the electrons have 
energies below 1 GeV, where their spectrum is steeper than 
N e o c p~ 2 , in order to obtain such a steep syn chrotron spectrum 
(see berezhko et all2003bHBerezhkol2005L for details). In ad- 
dition, the form of the X-ray synchrotron spectrum, including 
its cutoff, must be consistent with this magnetic field strength. 
This determines the effective field strength inside the SNR to 
within about 30 percent. In all cases considered in the past the 
field turned out to be amplified by a factor bet ween 5 and 10 rel - 
ative to the expected upstream field strength dVolk et all2005l) . 

The second method is also empirical and makes use of 
the synchrotron morphology of the SNR at keV energies 
dBerezhko & Volld l2004al IVolk et al J 12003) . Interpreting the 
thickness of sharp filaments in hard X-rays - typically in the 
line-free energy region between 4 and 6 keV - as the result 
of synchrotron losses during and after the electron acceleration 
process, we obtain an independent measure of the downstream 
magnetic field strength. Within an accuracy of 20 to 30 per- 
cent this field strength agrees with that determined by the first 
method. 
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The present status of the radio synchrotron measurements 
for SNR RX J1713.7-3946 does not determine the radio spec- 
tral index, and therefore we have only limited restrictions on 
the field strength. Fortunately, however, the synchrotron mor- 
phology gives a rather defi nite measure of the field. The re- 
cent XMM observations of Hira ga et"al] J2Q05I) show a radial 
X-ray profile (see their Figs.l and 2) which we interpret to 
show unequivocally the filamentary structure behind the outer 
SNR shock. The exponential angular width of the filament 
corresponds to 5.1 bins, of width Aifr = 25.6" each (Hiraga, 
private comm unication), within the radial profile shown by 
Hiraga et al. (2005). Therefore the profile width corresponds 
to Atfr « 2.2', or to a spatial scale L = 2x 10 18 cm at a distance 
of 1 kpc. 

As it was already demonstrated for other young SNRs 
jBerezhko et alJ l2003al iBerezhko & Volkl l2004at IVolk et alJ 
1200.4 . the measured width of the projected radial profile of 
the nonthermal X-ray emission determines the internal (down- 
stream) magnetic field according to the expression 

B A = [3mlc 4 /(4er 2 Q l 2 2 )] l/ \ Vl+c5 2 - 6T 2I \ (7) 

where 6 1 = 0.12[c/(r v)][V s /(o-c)] 2 , l 2 ~ L/7 is the ra- 
dial width of the X-ray emissivity q v (e v , r), is the classi- 
cal electron radius and cr is the total shock compression ratio. 
Substituting into this expression the values L — 2 x 10 18 cm, 
cr = 6.3 (see below), V s = 1800 km/s, and v = 1.8 x 10 17 Hz 
(e y = 0.7 keV), we obtain B d as 65 fiG. 

Given the finite resolution used to analyse these XMM data, 
which leads to a wider profile than the actual profile (see Fig|5] 
below), this value for Bj is clearly a lower lim it to the true 
amplified field. However, the Hir aga et al J (12005) XMM profile 
has the great advantage that it certainly corresponds to the outer 
shock, wher eas for the much narro wer 20" Cha ndra pr ofile 
measu red bv llJchivama et al] J2003I) . as analysed in lVolk et alJ 
J2005I) . the doubt remained, if that profile was due to the outer 
shock or rather due to some other discontinuity in the SNR 
plasma. On the other hand, there is at present only one such 
profile available, and we are left with the question as to how 
representative this profile and the derived field value is. From 
general gas dynamics considerations we argue that the proba- 
bility is extremely low that the XMM profile corresponds to a 
singular point at the shock surface. We therefore conclude that 
a lower limit to the magnetic field strength at the outer shock is 
65 [tG. Extrapolating from the 20" scale of the Chandra profile, 
the lower field limit could be as large as 230/iG. Examinations 
of other XMM or Chandra radial profiles are expected to con- 
firm this conclusion. 

Below we shall derive the effective field strength as a re- 
sult of fitting the synchrotron spectrum, in particular its X-ray 
part. For other SNRs this fit and the filament sizes give con- 
sistent results for the effective field. It will be shown that the 
spectral fit can be achieved with values Bj = 126 jjG as well 
as Z?d = 250 pG. These values roughly lie in between the two 
extremes discussed above. Such values of Bj are significantly 
higher than typical ambient magnetic fields, also for the con- 
sidered bubble wall in the molecular cloud (the density at the 
shock is still a factor ~ 100 lower than in the cloud which im- 
plies a strongly decompressed field upstream of the SNR shock 



compared to the value in the cloud). It must be attributed to 
field amplification at the shock front due to the strong wave 
producti on by the acceleration of CR s far into the nonlinear 
regime dBell & Lucekl200lUBel]|l2004 . 

3. Model 

Hydrodynamically, a core collapse supernova (SN) explosion 
ejects a shell of matter with total energy E SD « 10 51 erg 
and mass M e j equal to a few solar masses. During an ini- 
tial period the shell material has a broad distribution in ve- 
locity v. The fastest part of this ejecta distribut ion can be de- 
scribe d by a pow er law dM^/dv oc v - k fe.g. lJones et alll98li 
IChevalierll982l) . The interaction of the ejecta with the interstel- 
lar medium (ISM) creates a strong shock there which heats the 
thermal gas and accelerates particles diffusively to a nonther- 
mal CR component of comparable energy density. For core col- 
lapse SN explosions a value k — 8 appears appropriate which 
we adopt in this paper. 

As pointed out in the Introduction, our theory for this pro- 
cess is based on a fully time-dependent, spherically symmet- 
ric solution of the CR transport equations, coupled nonlinearly 
with the gas dynamic equations for the thermal gas component. 
Since all relevant equations, initial and boundary conditions for 
this model have already been described in detail in these pa- 
pers, we do not present them here and only briefl y discuss the 
mos t important aspe cts below (reviewed also by IVolkl d2003l) 
and lBerezhkol J2005I) ). 

The coupling between the ionised thermal gas (plasma) 
and the energetic particles occurs primarily through magnetic 
field fluctuations carried by the plasma which scatter ener- 
getic particles in pitch angle and energize them, especially 
in shock waves. The plasma physics of the field fluctuations 
is not worked out yet in full detail. However the accelerat- 
ing CRs effectively excite magnetic fluctuations upstrea m of 
the outer SN shock in the form of Alfven waves (e.g. Bell 
1978; Blandford & Ost rikerlll978h as a result of the stream- 
ing instability. In quasilinear approximation the wave ampli- 
tudes SB grow to very high amplit udes 6B > B, w here B 
is the average field strength JMcKenzie & Volk1ll982l) . Since 
these fluctuations scatter CRs extremely strongly, the CR dif- 
fusion coefficient is assumed to be as small as the Bohm limit 
K(p) = K(mc)(p/mc), where K{mc) = mc 3 /(3eB), e and m are 
the particle charge and mass, p denotes the particle momen- 
tum, and c is the speed of light. 

If Bism is the pre-existing field in the surrounding medium, 
then the strong streaming instability would suggest that the in- 
stability growth is restricted by some nonlinear mechanism to 
the level SB ~ Bism, any further turbulent energy being dissi- 
pated into the thermal gas JVolkfcMcKenziell 98n. The early 
attempts to give a full nonlinear descripti on of the magnetic 

field evolution in a numerical simulation ( Lucek & B ell 2000; 

II 1 

Bell & Lucek 2001) concluded that a considerable amplifica- 
tion of the "average" magnetic field in the smooth shock pre- 
cursor - produced by the finite CR pressure - should occur. 
Broadly speaking, it was expected that a non-negligible frac- 
tion of the shock ram pressure poV 2 (po is the ambient gas den- 
sity at the current SN shock position) is converted into mag- 
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netic field energy. Subsequently. lBel3 (2004) argued that this 
amplification is the result of a nonresonant instability, giving 
rise to what we shall call the effective upstream field Bq > Bi$ M , 
on top of which the Alfven waves grow to amplitudes SB ~ Bq. 
The Bohm limit and the gas heating are then to be calculated 
with B = Bo- 

In our analyses of the synchrotron spectrum of SN 1006, 
Tycho's SNR, and Cas A, such magnetic field amplifications 
were indeed found; they can only be produced as a nonlinear 
effect by a very efficiently accelerated nuclear CR component. 
Its energy density, consistent with all existing data, is so high 
that it is able to strongly excite magnetohydrodynamic fluctua- 
tions, and thus to amplify the upstream magnetic field Bism to 
an effective field Bq > Bi$ M , at the same time to permit efficient 
CR scattering on all scales, reaching the Bohm limit. The same 
large effective magnetic field turns out to be required, within 
the errors, by the comparison of this self-consistent theory with 
the morphology of the observed X-ray synchrotron emission, in 
particular, its spatial fine structure. We shall therefore also here 
allow for the possibility of an amplified field. 

The number of suprathermal protons injected into the ac- 
celeration process is described by a dimensionless injection 
parameter 77 <K 1 which is a fixed fraction of the gas parti- 
cles entering the shock front. For simplicity it is assumed that 
the injected particles have a velocity four times higher than the 
postshock sound speed. We have argued before that ion injec- 
tion is quite efficient at the quasiparallel portions of the shock 
surface, where it is ch aracterised by values 77 = 10~ 4 to 10" 3 
(see IVolk etaljEooi for details). Since this injection is ex- 
pected to be strongly suppressed at the quasiperpendicular part 
of the shock, one should renormalize the results for the nu- 
cleonic spectrum, calculated within the spherically symmetric 
model. The lack of symmetry in the actual SNR can be approx- 
imately taken into account by a renormalization factor / re < 1, 
roughly / re = 0.15 to 0.25, which diminishes the nucleonic CR 
production efficiency, calculated in the spherical model, and all 
effects associated with it. 

We assume that also electrons are injected into the accel- 
eration process at the shock front. Formally their injection mo- 
mentum is taken to be the same as that of the protons. Since the 
details of the electron injection process are poorly known, we 
chose the electron injection rate such that the electron:proton 
ratio K ep (which we define as the ratio of their distribution func- 
tions at all rigidities where the protons are already relativistic 
and the electrons have not been yet cooled radiatively) is a con- 
stant to be determined from the synchrotron observations. It is 
demonstrated below, that in the case of SNR RX J 17 13.7-3946 
the scarce existing data nevertheless give a possibility to esti- 
mate the values of Bq, r\ and K ep . 

The electron dynamics is exactly the same as that for pro- 
tons for electron rigidities corresponding to ultrarelativistic 
protons, as long as synchrotron losses are neglected. Therefore, 
beyond such rigidities and below the loss region the distribution 
function of accelerated electrons has the form f e {p) = K ep f{p) 
at any given time. The electron distribution function f e (p) de- 
viates only at sufficiently large momenta from this relation due 
to synchrotron losses, which are taken into account by supple- 



menting the ordinary diffusive transport equation by a radiative 
loss term. 

Clearly, from the point of view of injection/acceleration 
theory, we must treat K ep , together with Bo and 77, as a theoret- 
ically not very well constrained parameter to be quantitatively 
determined by comparison with synchrotron observations to 
the extent that they are available. 

The solution of the dynamic equations at each instant of 
time yields the CR spectrum and the spatial distributions of 
CRs and thermal gas. This allows the calculation of the ex- 
pected fluxes of nonthermal emission produced by the acceler- 
ated CRs. 

4. Results 

The computed gas dynamical characteristics of the SNR are 
shown in Fig^ They demonstrate in Fig^ that the calcula- 
tion fits the observed SNR size R s for the assumed distance and 
SNR age. The resulting shock velocity is V s « 1840 km/s. 

To obtain a good fit for the observed synchrotron and 
gamma-ray spectra (see below) we assume a proton injection 
rate 77 = 3 x 10~ 4 . We want to emphasise that the most reliable 
way to determine 77 is based on the measured spectral shape 
of the radio emission. Unfortunately this can not be done for 
the case of RX J1713. 7-3946 , because the shape of the radio 
spectrum is poorly known. The assumed injection rate is nev- 
ertheless in the same range as in other remnants, where this 
quantity can indeed be directly determined. It leads to a signif- 
icant nonlinear modification of the shock which, at the current 
age of t = 1612 yrs, has a total compression ratio <x = 6.3 
and a subshock compression ratio cr s - 3.3 (Fig^)- In the 
present case the shape of the TeV-spectrum (see below) is the 
only spectral characteristic which furnishes evidence that the 
SN shock is strongly modified. 

We note here that the lack of a restriction on the injection 
rate, resulting from the lack of knowledge of the radio spectral 
index, facilitates the fit of the spectrum. However, we will argue 
in section 5 that no basically different injection rate is possible. 

For its adopted density the hot wind bubble contains only 
a small amount of gas Mb ~ O.3M . This is consistent with 
the range of 0.07 < Mb < O.6M for the initial stellar masses 
M; in the range 15M < M ; < 2OM . Due to this fact the 
SN shock has deposited less than 20% of the explosion energy 
while propagating through the bubble (f < 700 yr), as seen 
from FiglQ;. Yet at the current epoch the SN shock has already 
swept up a considerable mass M sw m 2OM and therefore the 
ejecta has already transformed about 90% of its initial energy 
into thermal gas and CRs. 

In spherical symmetry the acceleration process is then char- 
acterised by a high efficiency: at the current time about 53% 
of the explosion energy have been transferred to CRs, and the 
CR energy content E c continues to increase to a maximum of 
about 60% in the later phase (FigUp), when the highest en- 
ergy particles start to leave the source. As typically predicted 
by the spherically symmetric model, such a CR acceleration 
efficiency is significantly higher than required for the average 
replenishment of the Galactic CRs by SNRs, corresponding to 
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Fig. 1. (a) Shock radius R s and shock speed V s ; (b) total shock 
(cr) and subshock (cr s ) compression ratios; (c) ejecta (E e j) and 
CR (E c ) energies as a function of time, normalised to the total 
mechanical SNR energy E sn . The vertical dotted line marks the 
current evolutionary epoch. 



E c m Q.lE sn . And as discussed above, the necessary devia- 
tion of the SNR from spherical symmetry requires a renormal- 
ization of the number of hadronic CRs. We choose the value 

/re = 1/5. 

With this renormalization the CRs inside RX J1713-3946 
SNR already contain 

E c = 0.53/ re £ sn * 2 x 10 50 erg. (8) 
The volume-integrated (or overall) CR spectrum 



N(p, t) = I6n 2 p 2 I drr z f(r, p, t) 



o 



(9) 




i°g(p/m p c) 

Fig. 2. Spatially integrated CR spectrum as function of particle 
momentum. Solid and dashed lines correspond to protons and 
electrons, respectively. 

has, for the case of protons, almost a pure power-law form 
N oc p~y over a wide momentum range from 10~ 2 /« p c up to the 
cutoff momentum p m . d!l » 2 x 10 5 m p c (Fig|2J. This value /? max 
is limited mainly by the finite size and speed of the shock, its 
deceleration and the adiabatic cooling effect in the downstream 
region (Berezhko 1996). As a result of the shock modification, 
the power-law index slowly varies from y = 2.2 at p <m v c to 
y = 1.7 at p 2* 10 3 m p c. 

The shape of the overall electron spectrum N e {p) deviates 
from that of the proton spectrum N(p) at high momenta p > 
p\ ~ 10 3 m p c, on account of the synchrotron losses during their 
residence time in the downstream region with a magnetic field 
strength m 130 pG which is assumed unifo rm in this region 
(B d = B 2 = o-Bo). Cf. lBerezhko et alJ J2002I) the synchrotron 
losses become important for electron momenta greater than 



P\ 

m D c 



1.3 



10 8 yr\/10yuG 
~B~T 



t 



(10) 



500 TOpC, in 



where p\ is relevant only in the evolutionary stage, when it be- 
comes lower than the electron cutoff momentum pf llax , p\ < 
Pmax- This is already the case at present. Substituting the SN 
age t — 1612 yr into this expression, we have p\ 
agreement with the numerical results (Fig|2|i. 

The shock continuously produces an electron spectrum f e oc 
p~ q , with q » 4, up to the maximum momentum p^ x (t) which 
is at the present time already much larger than p\. Therefore, 
within the momentum range p\ < p < p e mm , the electron 
spectrum is considerably steeper A^ e °c p~ 3 due to synchrotron 
losses taking place in the downstream region after the acceler- 
ation at the shock front. 

For p e m!a < p m ax the maximum electron momentum can be 
estimated by equating the synchr otron loss time an d the accel- 
eration time. This gives (e.g. lBerezhko et alJ2 002): 



= 6.7 x 10 4 



10 3 km/s 



(cr-1) (lOpG 



^cr(l+cr 2 )\ B 



(11) 
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Fig. 3. Spatially integ rated spectral energy d istribution of RX J1713.7-3946 . The AT CA radio data (cf . lAharonian et"ai1l2005l) . 
ASCA X-ray data (cf lAharonian et al.l2005l) . EGRET s pectrum of 3EG J 17 14- 3 857 JReimer & Pohl2002i) . CANGAROO data 
fenomoto et alJE()02l) . in red color) an d H.E.S.S. data jAharonian et alJl2005h . in blue color) are shown. The EGRET upper 
limit for the RX J1713. 7-3946 position JAharonian et alJl20 05) is shown as well (red colour). The solid curve at energies above 
10 7 eV corresponds to 7r°-decay y-ray emission, whereas the dashed and dash-dotted curves indicate the inverse Compton (IC) 
and Nonthermal Bremsstrahlung (NB) emissions, respectively. 



At the current epoch V s ~ 1840 km/s which leads to a maxi- 
mum electron momentum pf llax w 1.3 x 10 4 m p c, in agreement 
with the numerical results (Fig0. 

We note that during the last thousand years the SN shock 
speed has been going down rapidly. In previous evolutionary 
epochs it has therefore produced electron spectra with cutoff 
momenta /?^ ax larger than that of the current epoch. Due to this 
fact the overall electron spectrum has a relatively smooth cut- 
off, extending up to p ~ 10 5 m p c (see Fig|2}- Together with the 
synchrotron cooling it provides a very good fit of the observed 
X-ray spectrum (see below). 

The parameters K ep * 10~ 4 and Z?d = 126 /jG give 
good agreement between the calculated and the measured syn- 
chrotron emission in the radio to X-ray ranges (Fig0. The 
steepening of the electron spectrum at high energies > 10 3 m p c 2 
due to synchrotron losses and the smooth cutoff of the overall 
electron spectrum naturally yield a fit to the X-ray data with 
their soft spectrum. 

The overall broadband spectral energy distribution is 
displayed in Fig|3] together with the experimental data 
from ATCA at radio wavelength s , as estimated for the 
full remnant by lAharonian et alJ (2005), the X -ray data 
from ASCA, the EG RET spectrum o f the nearby source 
3 EG J1714-3857 JReimer & Pohllf2 002). and th e TeV y-ray 
spectra from CANGAR OO fenomoto et al J2002I) and H.E.S.S. 
( lAharonian et al.ll2.005h . The data also include the more recent 
EGRET upper limit based on the assumption that 3EG J1714- 



3857 is not physically associated with RX J1713. 7-3946 . The 
overall fit is impressive, noting that the choice of a few key 
parameters like 77, Z?d, and E sn in the theory allows a spectrum 
determination over more than 19 decades. The remainder of 
this section will give a detailed discussion of these spectra and 
of the morphology. 

In Fig|4] we separately present the synchrotron spectrum, 
produced at the current epoch by the accelerated electrons. For 
comparison we also include (through the dashed curve) a syn- 
chrotron spectrum, which would correspond to an artificial sce- 
nario with a proton injection rate so small (77 = 10" 5 ) that the 
accelerated nuclear CRs do not produce any significant shock 
modification and therefore also no magnetic field amplifica- 
tion. This corresponds to the test particle limit, and the low 
value Bo = 5 /jG was adopted for this case, for which syn- 
chrotron cooling is negligible. There are two meager differ- 
ences in the synchrotron spectra, corresponding to these two 
scenarios. The high-injection scenario leads to a steep radio 
spectrum S v oc v~ a with power law index a = 0.62 whereas in 
the test particle case a = 0.5. The quality of the existing radio 
data does not allow to distinguish these two scenarios. On the 
other hand the two spectra behave essentially different at X-ray 
frequencies v>10 1& . They demonstrate that only in the high- 
injection case case with its high, amplified magnetic field value 
fid ~ 100 pG the spectrum S v (v) has a smooth cutoff consis- 
tent with the experiment (see Fig@). In the test particle case the 



8 



E.G. Berezhko and H.J.Volk: Theory of CR production in SNR RX J1713.7-3946 



spectrum S v (v) has too sharp a cutoff to be consistent with the 
observations. 

The properties of small scal e structures o f SNR 
RX J1713.7-3946 seen in X-rays bv lUchivama et alJ d2003l) . 
and in particular bv lHiraea et al.l d2005l) . provide even stronger 
evidence that the magnetic field inside the SNR is indeed 
considerably amplified (see section 2.1). 

In order to find out if the filamentary structure found by 
fairaga et al.l2005h is indeed consistent with the high-injection 
case we show in Fig|3]the projected radial brightness profile 



J(e,p)oc J dxq(e,r — yjp 2 + x 2 ,x), (12) 

calculated for the X-ray energy e = e v = 1 keV. Here q(e, r) is 
the spectral luminosity of the nonthermal emission with photon 
energy e. The integration is performed along the line of sight 
x. For high injection the theory predicts the peak of the emis- 
sion just behind the shock front with thickness Ap/R s » 10~ 2 
that corresponds to an angular width Aif/ * 0.4'. This width 
is significantly thinner than the observed width of Atfr * 2.2'. 
However, one should take into account the smoothing proce- 
dure with which the data were obtained/presented. Instead of 
the actual profile J 7 (p) it gives a broadened profile 

dp'G( P ,p')J(p'), (13) 

oo 

through the (Gaussian) point spread function 

G(p,p') = (op V2M)exp[-(p-p') 2 /(2 C r2)], (14) 

where <x p = cr^d, and cr^, is the angular resolution of the in- 
strument (in radians). If we apply such a point spread func- 
tion with cr^ = 12.8" (corresponding to the smoothed XMM 
data) to the calculated profile, we obtain the broadened profile 
, which is shown in Fig|5]by the dashed curve. It has a width 
Ap/Rs w 6.6 x 10 2 , or At// w 2', and is consistent with the 
observational profile presented in Fig. 2 o f lHiraga et alJ (2005) 
which was discussed in section 2.1. 

In addition we present in Fig 13 the projected radial profile 
(dash-dotted line) which corresponds to the test particle limit. 
Since in this case the synchrotron losses are not relevant on ac- 
count of the much lower magnetic field value = 20 /iG, the 
peak of the emission is broader by a factor of five, inconsistent 
with the XMM observations. 

Morphologically the remnant has a rather complex struc- 
ture. This is particularly visible in the very bright western parts. 
Instead of a single, bright thin rim which should be identified 
with the actual position of the forward SNR s hock, there ar e 
two nested, narrow arc-like rims ( Cassam-Che nai et alJ l2004). 
The additional second inner bright rim can not be associated 
within our model with a structure situated inside the remnant at 
the corresponding distance behind the forward shock. On the 
other hand lCassam-Chenai et alJ J2004h found a positive corre- 
lation between the X-ray brightness and the absorbing column 
density. It is possibly explained by assuming that all bright 
spots observed in X-rays lie just behind the forward shock, 
whose shape is significantly distorted due to distortions of the 
dense, swept-up shell. Such shell distortions can in turn can 
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Fig. 4. Synchrotron photon flux density as a function of fre- 
quency. The solid line corresponds to the high-injection model, 
the dashed line corresponds to the test-particle approxim ation. 
The A TCA radio data and ASCA X-ray data ( Aharonian et alJ 
2005) are shown. 
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Fig. 5. Projected radial profile of the X-ray synchrotron emis- 
sion for the energy e v = 1 keV. The solid line corresponds 
to the high-injection model; the dashed line again represents 
the above profile, but smoot hed to the resolution of the XMM- 
Newton data used in Hiragaet al. (2005); the dash-dotted line 
corresponds to the test-particle limit. For purposes of presen- 
tation all profiles are normalised to their peak values in this 
figure. 

be expected because of inhomogeneities of the external ISM. 
In such a picture the outer rim, associated with the position of 
the forward SNR shock, is bright due to the projection effect, 
whereas the inner rim is bright because the corresponding part 
of the shock interacts with a denser medium. 

Concentrating now on y-ray energies the calculated non- 
thermal Bremsstrahlung, IC and 7r°-decay y-ray spectral energy 
distributions are shown in Fig |6] together with the existing ex- 
perimental data. (In this presentation the theoretical spectrum 
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Fig. 6. Nonthermal bremsstrahlung (dash-dotted line) IC 
(dashed line) and 7r°-decay (solid line) y-ray spectral en- 
ergy distribution as a function of y-ra y energy. The ob- 
served H.E.S.S JAharonian et alJ Eool and CANGAROO 
fenomoto et all2002h y-rav fluxes and t he EGRET spectrum of 
3EG J1714-3857 jReimer & Pohll2002l) and the EGRET upper 
limit are shown as well. 



can be dir ectly visually com pared with the H.E.S.S. results, cf. 
Fig. 20 of Ah aronian et al.l II2005I) .') According to the calcula- 
tion, the hadronic y-ray production exceeds the electron con- 
tribution by more than two orders of magnitude at all energies 
(Fig®. 

In the context of the CANGAROO observations of 
RX J17 13.7 -3946 , which covered only the northe rn part of 
the remnant fenomoto et all2002h . lButt et all J20Q2) as well as 
iReimer & Pohll J20 02) had argued that the GeV emission from 
the nearby unidentified EGRET source 3EG J1714-3857 is ei- 
ther associated with the SNR or an upper limit to the y-ray 
emission of RX J1713. 7-3946 . As one can see from Fig|6]such 
an upper limit does not contradict our calculated spectrum of 
RX J1713. 7-3946 . This is even true for an extrapolation of the 
observed H.E.S.S. spectrum in th e test-particle app roximation, 
with a spectral index of as 2 (I Aharonian et all2005l) . 

For energies e y < 3 TeV the theoretical y-ray spectrum is 
as hard as dF y /de y oc e y l& , whereas for e y > 10 TeV it has 
a smooth cutoff. We note that the y-ray cutoff energy e™ ax « 
0.1 cpmax is sensitive to the magnetic field strength B^, since 
the proton cuto ff momentum has a dependence p max oc /? s y s B d 
feerezhkoll 1 9961) . It is clearly seen from Fig|6]that the calcu- 
lated spectrum in a satisfactory way fits the H.E.S.S. measure- 
ments. 

Since the SN shock propagates through the shell with a ris- 
ing gas density, the TeV emission is expected to increase with a 
rate of about 0.4 %/yr, whereas the nonthermal X-ray emission 
goes down with a rate 0. 1 %/yr due to the shock deceleration. 

We note also that within our approach there is no need to 
suggest the interaction of the SN shock with the molecular 
cloud to be followed by a strong escape of high energy CRs, 
in order to explain the observed steepening of y-ra y spectrum 
at e r > 1 TeV, as proposed bv lMalkov et alJ J2005b . 

The projected radial y-ray emission profile, calculated for 
the energy e y — 1 TeV, is presented in Fig0 As a result of 
the large radial gradient of the gas density and the CR distribu- 
tion inside the SNR, the theoretically predicted radial profile of 
the TeV-emission is concentrated within a thin shell of width 
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Fig. 7. The y-ray emissivity for y-ray energies e y — 1 TeV 
as a function of angular radial distance. The calculated radial 
profile is given by the solid line; the dashed line represents the 
calculated profile smoothed with a Gaussian point spread func- 
tion of the width cr^ = 0.05°. For purposes of presentation both 
profiles are normalised to their peak values in this figure. 



Ap « 0.1R S . Due to the projection effect this width is seven 
times larger than the width of the 3-dimensional radial emissiv- 
ity profile, which is as thin as Ar m 0.01/? s . Since the H.E.S.S. 
instrument has a finite angular resolution we present in Fig0 
also the modified radial profile smoothed with the Gaussian 
point spread function with cr p = Ap = 0.08/? s , that corresponds 
to the angular resolution 2cr^, = 0.1°. As shown by Fig0the 
smoothed radial profile of the TeV-emission is much broader 
and - what is most interesting - it is characterised by a maxi- 
mum to minimum intensity ratio 7™ ax /7™ n = 2.2. Such a ratio 
is con sistent with the H.E.S.S. measurement ll Aharonian et alJ 
2005) which obviously gives only a lower limit to the sharpness 
of the y-ray profile. We conclude that the broad radial profile of 
the TeV-emission with 7™ x /7™ n * 2 measured by the H.E.S.S. 
instrument is indirect evidence that the actual radial profile is 
significantly sharper, with a higher ratio 7™ ax /7^ ln > 2.2. Also, 
comparing with the radial X-ray profile in Fig|5] the similarity 
of the radial profiles of the X-ray, as measured by ASCA and 
the radial TeV y-ray profiles, as obtained by H.E.S.S., must be 
seen as a consequence of the finite angular resolutions of these 
instruments. In reality, the X-ray profiles should be narrower 
than the y-ray profiles due to the synchrotron losses. 



5. Discussion and summary 

The theoretical results presented here are fully consistent with 
a dominantly hadronic origin of the observed TeV y-rays . 
Beyond that they of course imply that the overwhelming con- 
tribution to the nonthermal pressure comes from accelerated 
nuclear particles. Despite the observational fact that the broad- 
band emission is not fully spherically symmetric, with the 
brightness increasing from southeast to northwest, the quasi- 
spherical character of the X-ray rims on the western half of the 
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SNR suggests the approximate validity of a spherically sym- 
metric model, with clear limitations due to the non-spherical 
nature of the magnetic field configuration which are to be cor- 
rected for. 

The true difficulty for the theoretical description is there- 
fore not the complex geometry of the distribution of molecu- 
lar gas and X-ray absorbing column density, but rather the fact 
that several key parameters of this source are either not known 
or poorly constrained. This already concerns the distance and 
age of the object. We follow present consensus which puts the 
distance at 1 kpc and the age to about 1600 years. We also 
concur with other authors and argue that the primary explosion 
must have been a type Il/Ib SN event with a massive progen- 
itor star whose mass loss in the main sequence phase created 
a hot wind bubble in a high-density environment. However, at 
present these remain assumptions, as well-founded as they are. 
The solution for the overall remnant dynamics then yields the 
value for the expansion velocity of the outer shock, given the 
total mechanical energy E sn released in the explosion. To obtain 
a consistent solution for the broadband nonthermal emission 
we chose E sn = 1 .8 x 10 51 erg which is already on the high side, 
also for core collapse SNRs with initia l main seque nce progen- 
itor masses between 15 and 20M Q llHamiivl l2003). Therefore 
this particular parameter value is limited from above. 

The most difficult aspect is the absence of a reliable radio 
synchrotron spectrum. In other SNRs, such as SN 1006, Cas A 
and Tycho's supernova, the deviation of the radio photon index 
a from its test particle value a = 0.5 allows the determination 
of the injection rate of the dominant nuclear particles into the 
acceleration mechanism at the forward SNR shock, and thereby 
the nonlinear modification of the shock and of the acceleration 
process. Fitting then the overall synchrotron spectrum, from 
radio to X-ray frequencies, at the same time the effective mag- 
netic field inside the SNR can be determined with reasonable 
accuracy. 

In the absence of a radio spectral index we have estimated 
the strength of the effective field by analysing the observed 
width of an X-ray filament in a smoothed radial profile which 
clearly indicates the position of the outer shock. The result is 
a lower limit to the effective magnetic field strength of about 
70 /jG. Taking a value of about 130 /iG to be representative 
for the postshock value over 80% of the remnant's outer sur- 
face, we can calculate the IC y-ray emission from the known 
X-ray synchrotron emission. Its negligible contribution - at 
the present epoch - to the observed TeV y-ray flux, simi- 
lar to that from nonthermal Bremsstrahlung, shows that SNR 
RX J1713. 7-3946 must be a hadronic y-ray source. This is our 
first major conclusion. 

The above magnetic field strength is clearly amplified with 
respect to its expected value at the inner border of the hot wind 
bubble. Theoretically this amplification can only come from 
an efficiently accelerated nuclear particle component. And in- 
deed, choosing as a consequence a proton injection rate sim- 
ilar to those determined for the other young SNRs mentioned 
above, it is possible to consistently describe the entire nonther- 
mal broadband spectrum over 19 orders of magnitude in photon 
energy, including now the observed 7r°-decay y-ray spectrum. 
This is our second important conclusion, and it is far from triv- 
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Fig. 8. Nonthermal Bremsstrahlung (dash-dotted line), IC 
(dashed line) and 7r°-decay (solid line) y-ray spectral energy 
distributions as a function of y-ray energy for the case, where 
Eq^]is imposed as an additional constraint. The following pa- 
rameters have been modified compared to those on which Fig|6] 
is based: E sn = 2.5 x 10 51 erg, M ej = 5M G , B = 50 yt/G, 
rj = 10~ 4 , which give rise, at the current epoch, to the values: 
V s = 1651 km/s, cr = 4.93, and cr s = 3.71. N g was chosen as 
N g = 1.76 crrT 3 , in order to be able to visually compare the 
form of the resulting hadronic y-ray energy spectrum with the 
observations. The data plotted are the same as in Fig[6] 

ial. On the other hand, one could argue that this remarkable re- 
sult could be achieved because the observational constraints are 
not as numerous as they would be in the presence of a complete 
multi-wave-length knowledge. It therefore remains to justify 
the high nuclear injection rate to positively predict the ampli- 
tude of the y-ray flux. A reasonably good radio spectrum could 
achieve this. The possibility of measuring a good thermal X-ray 
spectrum seems very limited from the known XMM-Newton 
results. 

In the following we shall give a further going argument 
which is based on the physical idea that the magnetic field am- 
plification is driven by the gradient of the CR pressure upstream 
of the outer shock. In fact, for all the thoroughly studied young 
SNRs mentioned above, the ratio of magnetic field energy den- 
sity Bq/87t in the upstream region of the shock precursor to the 
CR pressure P c is about the same. Within an error of about 50 
percent we have 

bI/(SttP c ) * 5 x 10~ 3 . (15) 

Requiring now this semi-empirical relation to hold also for 
SNR RX J1713. 7-3946 , we can ask ourselves, what magnetic 
field strength, larger than the lower limit which we could de- 
termine from the XMM observations, can reproduce the ob- 
servations. For this consideration we assume in addition that 
the total mechanical explosion energy E SB is not larger than 
2.5x 10 51 erg, i.e. that RX J 17 13.7-3946 is not an outlier in the 
distribution of core collapse explosion energies with main se- 
quenc e progenitor masses below 2OM (compare also lHamuvl 
120031) . Together with other parameters (see Fig|8}, roughly 
equal to those given in section 2, even this large value of 
E sn can reproduce the observations, not counting here the y- 
ray observations. The associated proton injection rate is then 
also determined. Finally, the resulting theoretical hadronic y- 
ray spectrum must be compared to the TeV observations. 
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Obviously the above additional constraints make it more diffi- 
cult to fulfil the sum total of the constraints. In addition, the rel- 
evant gas density at the bubble wall remains poorly constrained. 
However, the ratio E c /E sn of nuclear CR energy to mechanical 
explosion energy is now limited from below. 

In a limited coverage of the entire parameter space it turns 
out that (i) Bd * 250 fiG (ii) that the ion injection rate rj = 1CT 4 
is a factor of 3 smaller than assumed in the previous sec- 
tion, and (iii) that the total CR energy inside the remnant is 
E c = 0. 22 f ie E SB as 10 5() erg. It will grow in the spherically 
symmetric case to about 50% of E sn in another few thou- 
sand yrs, in general agreemen t with earlier studies of wind- 
SNe (Berezhko & Volk 2000a). Choosing a present-day pre- 
shock density which is about 50 percent higher than assumed 
in section 2, one can even fine-tune the amplitude factor of 
the hadronic y-ray flux to the observed amplitude, in order to 
compare with the observed spectral form. The result is shown 
in FigfJ] It demonstrates that the new solution is less nonlin- 
early modified than before. However, the observed part of the 
spectrum is quite well reproduced. It is quite clear from Figs|6| 
and|8]that measurements of the y-ray spectrum at low energies 
e y S\0 GeV, which can be done with the y-ray space telescope 
GLAST in the future, could provide very useful information 
about the degree of SNR shock modification. 

We conclude that the present observational knowledge of 
SNR RX J1713. 7-3946 can be interpreted by a source which 
ultimately converts about 10% of the mechanical explosion en- 
ergy into nuclear CRs and that the observed high energy y-ray 
emission of SNR RX J 17 13.7-3946 is of hadronic origin. 
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